99 wetlands

Retrieving wetland boundaries from the National Wetlands Inventory (NWI)

The National Wetlands Inventory (NWI) is a comprehensive geospatial inventory of U.S. wetlands. It is a publicly available resource that provides detailed information on the abundance, characteristics, and distribution of wetlands in the United States. The NWI dataset is maintained by the U.S. Fish and Wildlife Service (USFWS) and is used by a wide range of stakeholders, including federal, state, and local agencies, researchers, and conservation organizations.

The notebook demonstrates how to retrieve wetland boundaries from the NWI using the leafmap Python package. The NWI dataset is available as a web service, which allows users to access the data programmatically. The leafmap package provides a simple interface for querying the NWI web service and visualizing the wetland boundaries on an interactive map.

Uncomment the following line to install leafmap if needed.

In [1]:

# %pip install -U leafmap

# %pip install -U leafmap

In [2]:

import leafmap

import leafmap

In [3]:

m = leafmap.Map(center=[47.223940, -99.885386], zoom=14)
m.add_basemap("HYBRID")
m

m = leafmap.Map(center=[47.223940, -99.885386], zoom=14) m.add_basemap("HYBRID") m

Google Maps API key is required to use Google Maps. You can generate one from https://bit.ly/3sw0THG and use geemap.set_api_key(), defaulting to Esri basemaps.

Out[3]:

Using point geometry¶

In [4]:

point_geometry = {"x": -99.907986, "y": 47.216359}

point_geometry = {"x": -99.907986, "y": 47.216359}

In [5]:

gdf = leafmap.get_nwi(point_geometry)

gdf = leafmap.get_nwi(point_geometry)

In [6]:

m = leafmap.Map(center=[47.223940, -99.885386], zoom=14)
m.add_basemap("HYBRID")
m.add_nwi(gdf, add_legend=True)
m

m = leafmap.Map(center=[47.223940, -99.885386], zoom=14) m.add_basemap("HYBRID") m.add_nwi(gdf, add_legend=True) m

Google Maps API key is required to use Google Maps. You can generate one from https://bit.ly/3sw0THG and use geemap.set_api_key(), defaulting to Esri basemaps.

Out[6]:

Using bounding box¶

In [7]:

bbox_geometry = {"xmin": -99.9023, "ymin": 47.211, "xmax": -99.8556, "ymax": 47.2325}

bbox_geometry = {"xmin": -99.9023, "ymin": 47.211, "xmax": -99.8556, "ymax": 47.2325}

In [8]:

gdf = leafmap.get_nwi(bbox_geometry)

gdf = leafmap.get_nwi(bbox_geometry)

In [9]:

m = leafmap.Map()
m.add_basemap("HYBRID")
m.add_nwi(gdf, layer_name="Wetlands", zoom_to_layer=True)
m

m = leafmap.Map() m.add_basemap("HYBRID") m.add_nwi(gdf, layer_name="Wetlands", zoom_to_layer=True) m

Google Maps API key is required to use Google Maps. You can generate one from https://bit.ly/3sw0THG and use geemap.set_api_key(), defaulting to Esri basemaps.

Out[9]:

Using GeoDataFrame¶

In [10]:

bbox = [-99.8653, 47.3952, -99.7498, 47.4441]
bbox_geometry = leafmap.bbox_to_gdf(bbox)
bbox_geometry

bbox = [-99.8653, 47.3952, -99.7498, 47.4441] bbox_geometry = leafmap.bbox_to_gdf(bbox) bbox_geometry

Out[10]:

	geometry
0	POLYGON ((-99.8653 47.3952, -99.8653 47.4441, ...

In [11]:

gdf = leafmap.get_nwi(bbox_geometry)

gdf = leafmap.get_nwi(bbox_geometry)

In [12]:

m = leafmap.Map()
m.add_basemap("HYBRID")
m.add_nwi(gdf, layer_name="Wetlands", zoom_to_layer=True)
m

m = leafmap.Map() m.add_basemap("HYBRID") m.add_nwi(gdf, layer_name="Wetlands", zoom_to_layer=True) m

Google Maps API key is required to use Google Maps. You can generate one from https://bit.ly/3sw0THG and use geemap.set_api_key(), defaulting to Esri basemaps.

Out[12]:

Using HUC8 watershed boundary¶

In [13]:

huc8 = "11120104"

huc8 = "11120104"

In [14]:

gdf = leafmap.get_nwi_by_huc8(huc8, quiet=False)

gdf = leafmap.get_nwi_by_huc8(huc8, quiet=False)

Downloading...
From: https://documentst.ecosphere.fws.gov/wetlands/downloads/watershed/HU8_11120104_Watershed.zip
To: /tmp/HU8_11120104_Watershed.zip

  0%|          | 0.00/4.36M [00:00<?, ?B/s]

 60%|██████    | 2.62M/4.36M [00:00<00:00, 8.67MB/s]

100%|██████████| 4.36M/4.36M [00:00<00:00, 14.0MB/s]

Extracting files...

In [15]:

m = leafmap.Map()
m.add_basemap("HYBRID")
m.add_nwi(gdf, layer_name="Wetlands", zoom_to_layer=True)
m

m = leafmap.Map() m.add_basemap("HYBRID") m.add_nwi(gdf, layer_name="Wetlands", zoom_to_layer=True) m

Google Maps API key is required to use Google Maps. You can generate one from https://bit.ly/3sw0THG and use geemap.set_api_key(), defaulting to Esri basemaps.

Out[15]: